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Using a soft pion theorem based on chiral symmetry and a A (1232) resonance 
model we propose an estimate for the production cross section of low energy pions 
in the deeply virtual Compton scattering (DVCS) process. In particular, we express 
the ep — ► e'y-KN processes in terms of generalized parton distributions. We provide 
estimates of the contamination of the ep — > cyp DVCS observables due to this asso- 
ciated pion production processes when the experimental data are not fully exclusive, 
for a set of kinematical conditions representative of present or planned experiments 
at JLab, HERMES and COMPASS. 
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I. INTRODUCTION 



The past few years have witnessed an intense theoretical activity in the field of Generalized 
Parton Distributions (GPDs). They parametrize thenon perturbative content of the hadrons 
and are essentially matrix elements of the form ^ , 



(N'\q a (x») q/3 (0)\N), (1) 

where x M is a light-like vector {x 2 = 0). In (JTJ q a is the quark field with Dirac index a and 
N, N' are two hadronic states which can differ either by their structure or by the kinematics. 
The limit N = N' gives the ordinary parton distributions. 

In the Bjorken limit the GPDs appear in the amplitudes of exclusive reactions such as 



f*(q) + N — M + N', 



(2) 



where 7* is a highly virtual photon (Q 2 = —%cf — > 00 ) and M can be a meson or a photon, 
in which case the reaction is called Deeply Virtual Compton Scattering (DVCS). It is widely 
believed that DVCS is conceptually the cleanest process to access the GPDs. The basic 
reason is that the amplitude for meson production involves the not so well known meson 
wave function. By contrast the final real photon in DVCS can be considered as ppintlike 
because photon structure effects (VDM-like) are suppressed by powers of l/Q 2 

This theoretical simplicity of DVCS is to some extent counterbalanced by its experimental 
difficulty which is certainly greater than in the case of meson production. Nevertheless a 
number of experimental attempts to measure photon electroproduction off the proton have 



been performed or are in progress 



For the interpretation of these experiments 



to be fruitful it is compulsory to have a control on the exclusivity of the final state. If we 
consider DVCS on the proton, the most interesting final state is the proton itself, 

7* + P -> 7 + P (3) 

and from now on we restrict our attention to this case. For further reference we call it 
"Elastic DVCS". Q 

The problem for most experiments is that the experimental energy resolution is not 
good enough to isolate the exclusive (7 + p) channel because it is separated from the first 
strong inelastic channel by only the small pion mass. In practice such data will always be 
contaminated by the reactions, referred to as Associated DVCS (ADVCS) : 

7*+p-^7 + p + 7T° or 7* + p — > 7 + n + n + , (4) 

where 7r +, ° is a low energy pion which escapes detection. Note that experiments which have 
sufficient energy resolution to distinguish the 'jnN final state from the jp one can study the 
process (JH) for itself [ljj, thereby enlarging the scope of virtual Compton scattering. This 
is likely to be the case of the experiments planned at Jlab. 

In this paper we propose a calculation of the cross section for reaction (jlj) and we com- 
pare it with the elastic channel (JHJ) by integrating over the pion momentum up to a given 
cutoff. As the dangerous pions are those which have a small energy we use the soft pion 
techniques based on chiral symmetry. In the soft pion limit, that is k n — > where k n is the 



pion 4-momentum, this allows to evaluate the associated DVCS amplitude using the same 
generalized parton distribution as in the elastic case. So, in a relative sense, this evaluation 
is model independent. The inherent limitation of this approach is that, apart the chiral 
limit (m n — > ), it gives a reliable estimate only for a small center of mass (CM) energy of 
the final pion-nucleon pair. Typically the upper limit is set by the excitation energy of the 
first A resonance, that is about 300 MeV. In order to increase the range of validity of our 
estimate we propose a (model dependent) estimate of the associated DVCS corresponding 
to 

7 * -\-p — > 7 + A. 

For this we use the large N c limit as a guidance which allows one to relate the GPDs of the 
N — > A transition to the ones of the N — > N transition 

In practice one measures the (I, I', 7) reaction where / denotes either a muon or an 
electron. The amplitude T 11 '' 1 for this reaction is the coherent sum of T vcs , the virtual 
Compton scattering amplitude and of T BH the amplitude of the Bethe-Heitler process where 
the final photon is emitted by the lepton. When the reaction produces only a proton (elastic 
case) the calculation of this "elastic-BH" amplitude only involves the elastic form factors 
of the proton, which are well known. When a pion is produced together with the final 
nucleon, then the corresponding "associated-BH" (ABH) amplitude involves the pion electro- 
production amplitude. For consistency reasons we shall evaluate this amplitude in the same 
framework as the associated DVCS. 

Our paper is organized as follows: in Section [n] we specify the kinematics and give 
the relevant expressions for the amplitudes and cross sections. In Section ITTT1 we remind the 
leading order approximation for the DVCS amplitude and define the twist-2 quark operators 
which control the DVCS amplitude both in the elastic and associated case. In Section Hvl we 
remind the definition of the GPDs in the elastic case. In Section El we derive the soft pion 
theorems relevant for the evaluation of the associated DVCS and BH amplitudes. In Section 
IVII we present a model to estimate the associated DVCS and BH in the A region. Section 
IVHI is devoted to a presentation and discussion of our results in kinematical situations of 
interest. Section lYUH is our conclusion. 
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Figure 1: Graph of the virtual Compton process 

II. PRELIMINARY 

In Table |U we have collected the characteristics of the particles involved in the reaction. 
We define q as the momentum of the virtual photon exchanged in the VCS process, that is 
q = k — k! . So the BH virtual photon has q — q' . The relevant Lorentz scalars are 



Q = -q , t y = (q-q) , x B 



2p.q ' 



W 2 = (p + q - q 



l\2 



(5) 



We shall also note t — (p' — p) 2 which coincides with £ 7 in the elastic case or in the limit 

K = o. 

In the one photon exchange approximation and in the Lorentz gauge we have the generic 
expressions for the VCS and BH amplitudes : 



rpVCS 

jiBH 



±e*e';H^—u{l') lv u{l), 



(6) 



-e 3 J v - u(l') 
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Figure 2: Graphs of the Bethe-Heitler process 
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7.(fc - q') - mi 



u{l) 



T 



ZZ'7 



VCS | rpBH 



(7) 
(8) 



They correspond to the graphs of FigJT] and El respectively. In the above expressions the 



charge e = y47r/137 has been factored out and the ± sign refers to the charge of the lepton 
beam. The spinors are normalized to uu = 2m with m a generic mass and the photon 
polarisation is normalized to e'.e' = —1. The amplitudes H and J are defined by 



= -i< f\J d 4 xe- iq - x T[j"(0),f(x)} \i >, 
J" = </|j"(0)|i>, 



(9) 
(10) 



where i, f are the appropriate hadronic states and j(x) is the electromagnetic current carried 
by the quarks 

f( x ) = qYQq(x) = J2 QfQfYqAx), (ii) 

/ 

with Q the diagonal charge matrix Q= [2/3, — 1/3, • • •]. The matrix element of the electro- 
magnetic current between nucleon states has the usual form factor decomposition : 



< N(p')\f(Q)\N(p) >=u(p') 



Ff' n (t)-f + iF£' n (t) 



2M 



u(p) 



(12) 



where p, n refer to the proton or the neutron and t = (p — p') 2 - With our normalizations we 
have Ff (0) = 1, Ff(0) = 0, (0) = 1.79, F™(0) = -1.91. 



We shall note T e '[ 7 the amplitude for a final hadronic state containing only a proton and 
^os 7 ( 7r+ '°) the amplitude for producing a nucleon and a pion of charge +1 or 0. A similar 
notation will be used for the cross sections. 

To fix the idea, we consider an experiment where (k',q') the momenta of the final lepton 
and photon are measured while the final hadronic state p' or (j/, k n ) is not observed. The 
elastic (ll'j) event is selected by the missing mass condition 

W 2 = (k - k' + p - q'f = M 2 . (13) 



In pratice the events will be integrated up to a cutoff W max which generally exceeds the pion 
production threshold. In the following we assume that the experimental resolution on W is 
nevertheless good enough so that the production of more than one pion can be neglected. 
If we assume that the one particle states are normalized as : 



< k\k' >= (2ir) 3 2k 6(k-k'), 



(14) 



we have the following expressions for the invariant cross sections : 



da"' 7 = d 4 T5(W - M)dW 



kl 



u'j 2 



el. 



doZ!(* + >°) = d^dW^ J dk; \T^{^) 
In the above equations the common phase space factorj^] is 



(15) 
(16) 



d 4 r 



dxsdQ 2 dt 1 d& 



\2%{2Tty{j>.k) 2 x B Jl + 4M 2 x 2 B /Q 2 



(17) 



where $ is the angle between the planes (k, k') and (q,q f ) and k* = k^ k* stands for the 
pion momentum in the frame defined by p + q — q' = 0, that is the rest frame of the final 
pion-nucleon pair. One has 



kl 



2 _ W 4 - 2W 2 {M 2 + ml) + (M 2 - ml) 2 
~ W 2 ' 



(18) 



From Eqs. (|15fl6j) we see that the contamination of the elastic reaction by the associated 



one is measured by the dimensionless quantity 



K(Q 2 ,x B ,t y ,$) 



dW 



I dk* J2hel. 



rpll'-yi 



71 



+ 



167T 3 



I let 



P 



(19) 



where J2hei denotes the sum over all the helicities. This is the quantity we want to evaluate. 



III. HANDBAG APPROXIMATION FOR THE VCS AMPLITUDE 

We need to evaluate the hadronic tensor (Eq|H|) in the generalized Bjorken limit 
which here is defined by : 

Q 2 -> oo, x B fixed, t 7 /Q 2 -> 0, W 2 /Q 2 -> 0, (20) 

where the two last conditions are necessary to have factorisation of the amplitude The 
condition that W remains small with respect to Q is not explicitely stated for elastic VCS 
since W = M in this case, but it is necessary for the associated VCS. The factorisation 
theorem garantees that the amplitude factorizes in a hard part which can be computed in 
perturbation theory as a series in as(Q 2 ) and a soft part which depends on the long distance 
structure of the hadron. In this work we consider only the leading term of the hard part 
which amounts to evaluate the amplitude in the handbag approximation as shown in Fig. El 

For this purpose we adapt the formulation of [l| to our problem. In the Bjorken limit 
both the elastic and associated VCS are light-cone dominated. Therefore it is convenient to 
introduce two light-like vectors n and p to record the flow of hard momentum. These Sudakov 
vectors are chosen to be in the hyper-plane defined by the virtual photon momentum q and 
another vector P which is related to the target-ejectile motion. In our case it is convenient 
to choose : 

P = for the elastic VCS, (21) 

P = P + P 2 +fc " for the associated VCS, (22) 



Figure 3: The handbag approximation to the VCS amplitude. The crossed term is not shown. 



so that, for both reactions, P has the same expression in terms of p, q, q' that is : 



(23) 



Other choices of P are allowed, for instance P = p, but the formulation is more symmetric 
if one uses the definition (|23|) . One imposes the normalization conditions : j^] 



n 2 = p 2 = 0, n.p = 1, n.P = 1, 



(24) 



and one defines the variable £' as : 



Q 2 P.q 



AP.q 2P 2 



-1 + 



1 + 



Q 2 P 2 



\ (P-QY 



x b /2 
1 - x B /2 



in the Bjorken limit. (25) 



From this one gets the decomposition : 



P 

Q 



p2 

2 



p + —n, 



Q 2 



(26) 



where 



p2 = w M!±Z!_i 

2 4 



(27) 



An arbitrary vector a has the covariant decomposition 



a = a.np + a.pn + a(_L), (28) 

with the _L component defined by a(_L).n = a(_L).p = 0. If one restricts to the case where 
the final photon is real one then gets : 

/ W 2 - M 2 \ 
A = -2£p+ £P 2 + U + A(±), 



where 



^ Q 2 -t,-2?(W 2 -M 2 ) 
? = ? g 2 + 4 ^/2p2 ^ 6 m the Bjorken limit. 



If one keeps the leading light-cone singularity of the time ordered product of currents in 
Eq. (jnj one has, at leading order in a s {Q 2 ) (3| : 

H^= fdxC!?(x,t)<f\52& / Sf(x,n)\i> 

f 

+ J dxC^(x,0 < f\^Q 2 S 5 f (x,n)\i > +0 (j^j , (29) 

where the hard scattering coefficients (in which we consistently make the approximation 
£ ps £') write : 

<T = + n"^ 1 - <T) f ^T- + _J . ) , (30) 

2 yx — t, + ie x + ^ — isJ 

CT(x,0 = f ^— —7 — -1 , (31) 

2 H \x — £ + ie x + £ — is J 

The twist-2 operators Sf and Sj, for which we often use the global notation or are 
defined by : 

r d\ 

S f (x,N) = J —e iXx q f (-\N/2)"/.Nq f (\N/2), (32) 
r d\ 

Sj(x,N) = J —e iXx q f (-XN/2)j.Nj h q f (XN/2) (33) 



P 2 -Pl 



P2-P1 




Figure 4: Representation of the twist- 2 operator. 



with N an arbitrary light-like vector. They obviously satisfy the scaling law : 



(34) 



In Eq. (12911 one has N — n, which is the only relic of the hard scattering in the soft matrix 
elements. 

In Fig. 131 we have represented the matrix element of between generic hadronic states 
i(pi) and ${3)2)- The quark lines are on mass shell Fock states quantized at equal light cone 
time x^N^ 1 = 0. If one labels the momenta of the initial and final active quarks as K — (p 2 — 
Pi)/2 and K + (p 2 — Pi)/2 respectively, then the integration over A in Eqs. (|32"33|) implies 
that x = K.N. On the other hand, from the additivity of the longitudinal momentum, one 
has : 



(p 1 +p 2 ).N 



< K.N < 



(Pl+P2)-N 



2 - ' - 2 

So, if one chooses the normalization N.(pi + p 2 )/2 = 1 as in (12*41) . the integration over x is 
effectively restricted to the interval [—1,1] in Eq. "2EJ). 



In Eq. the initial state is always the nucleon, with momentum P 



p. On the 



other hand the final state, which has momentum P + -j = p + A, can be either the nucleon 
(p + A = p') or the pion-nucleon system (p + A = p' + k n ). Note that for given p, q, q' 
the Sudakov vectors are the same in both cases. This is not essential but it simplifies the 
formulation. 

Finally, we point out that because expression (jSHj) is accurate up to terms of order 1/Q 
gauge invariance is not strictly satisfied. One finds instead : 



g^~A(i.), H^q„ = 0. 



However it has been shown in [1J] that electromagnetic gauge invariance is restored by adding 
to a term which is explicitely of order ^J— t 7 /Q, and therefore does not contribute in 
the Bjorken limit. The effect of this extra term has been calculated for realistic situations 
and found to be negligible So gauge invariance is not a serious issue as long as the 
conditions ((25}) are satisfied. 



IV. GENERALIZED PARTONS DISTRIBUTIONS FOR THE ELASTIC CASE 



The matrix elements of are the non perturbative inputs of the calculation. If we 
note N(pi),Njj) 2 ) two generic nucleon states, the elastic matrix elements are parametrized 
according to 



-VIP2, 
Q: 

<N(p 2 )\S f {x,N)\N(p 1 )>-- 



1 r; u (P2) 



P 12 .N 



Hf/Nl-N + iE f/N 



2M 

<N{p 2 )\S 5 f (K,N)\N(pi) >= 

N.(p 2 -Pi) 



u(pi), (35) 



1 rMP^ 



PiaJV 



H f / N j.N + E f/N - 



2M 



l 5 u(Pi), (36) 



where we note P\ 2 = (pi+j>2)/2. The Generalized Parton Distributions (GPDs) H, E, H, E 
are a priori functions of x and of the invariants which can be formed with p±, p 2 and N, that 
is : 

GPD = GPD {x, (p 2 - Pl ).N, P 12 .N, t 12 = ( Pl -p 2 ) 2 } . 
However from the scaling law (f3i|) and the definitions (|35|36|) we have : 

GPD {lux, u(p 2 - Pl ).N, uP 12 .N, t l2 ] = GPD {x, (p 2 - Pl ).N, P 12 .N, t 12 ] 

so the GPDs are only function of t\ 2 and of the ratios x/P\ 2 .N and (p 2 —pi).N/Pi 2 .N . The 
standard choice for the functional dependence is : 

GPD{x, (p 2 - Pl ).N,P 12 .N, t 12 ) 

■°™{— ^• € - J ^S £ '-}- (37) 



Of course in the elastic case one has x = x i2 and £ 12 = £ but when one of the nucleons is an 
intermediate state which is off the energy shell, as is happens in the ADVCS process (see 
subsection IV Djl . this is no longer true. 

The interest and properties of the GPDs are presented in several reviews (lfil. Ilfil IttI. . 
In Section IVHI we explain the model used for our estimates. Note that one can devise a 
parametrization analogous to (|35|36|) when the final state is the ttN system [3], but for 
practical purposes that is not very useful as it involves many new unknown functions. 

For further use we introduce iso-vector GPDs through : 

S v , a (x,N) = I ^e lA ^(-AiV/2)7.iV^n ud g(AiV/2), (38) 

where (r a , a = 1, 2, 3) are the Pauli matrices and 11^^ project on the u, d part of the quark 
flavor multiplet. The parametrisation in term of GPDs then writes 

< N(p 2 )\S v , a (x,N)\N( Pl ) >= u(p 2 ) [H vT N+- ■ ■} -fu( Pl ) (39) 

where now the Dirac spinor is also a iso-spinor for the nucleon and t n are the Pauli matrices 
acting on it. A similar definition holds for Sy. The fact that (H v , E v , ■ ■ .) are independent 
of a and of the isospin of the initial and final nucleons is a consequence of isospin symmetry. 
As a consequence of their definitions, we have the following (generic) relations : 

Hy = H u / P — H d / P = H d / n — H u / n . (40) 



The iso-vector GPDs are thus controlled by the valence quarks in so far as the sea is flavor 
symmetric. The parametrisation defined by Eqs. (I38I39|) is necessary to describe reactions 
where the nucleon charge is changed, as is the case when a charged pion is emitted. 



V. SOFT PION THEOREMS 



A. Soft pion expansion 



To evaluate the amplitude for the reaction : 



l + N-^l' + >y + N + -K, 



we see from Eqs (|6|7ll()|29j) that we need the matrix elements : 



<iV7r|i At |A^>, 



(41) 



which is the pion electro-production amplitude, and 



< Nir\Y l &fSf ) \N>, 



(42) 



/ 



which we call (loosely) the pion twist-2 production amplitudes. To simplify the notations 
we do not write the isospin label of the pion or of the various iso-vector quantities if it is 
not necessary. 

We first recall the main steps in the derivation of a soft pion theorem for a matrix element 
of the form : 



where A is either the electromagnetic curre nt . Eq. (|4T"j) . or one of the twist-2 operators, 
Eq. (j42l) . In the approach of de Alfaro et al. (l3,l20|, which we follow here, one considers a 
physical pion at rest and looks for the dominant contribution as m n — ► 0. In the chiral limit, 
the isovector axial current : [441 




< JV7r(0)|i4|JV > 




jl' a = 97^75 yh^g, a = 1, 2, 3 



(43) 



is conserved. Therefore the axial charge operator : 




(44) 



is time independent. Due to the spontaneous breaking of the SU(2) L x SU(2) R symmetry 
down to SU(2)v, the axial charge does not anihilate the vacuum state but instead creates 
states with zero momentum (massless) pions. This is the origin of the soft pion theorems. 
The symmetry is also explicitely broken by the small quark mass and for our purposes it is 
sufficient to use the Partially Conseved Axial Current (PCAC) formulation (2l| which we 
write in the form : 



To derive the low soft pion theorems one does not need to be more specific about the non 
conservation of the axial current . One just needs to know the transformation properties 
of the relevant operators under the symmetry group and that the symmetry breaking is of 
order m%. To exploit these hypotheses one first defines the iso- vector operator :|2G] 



«=i* + £s«'L/ <45 » 



which satisfies 



< nP(K)\Q a \0 > = 0, (46) 

< 0\Q a \^(K) > = 2 i f n E 7T (K)(2nf5(a,(3)5(K), (47) 



where E n (K) = y K 2 + and is the pion decay constant. A general matrix element of 
Q obviously satisfies : 

~™ 

{E 1 - E 2 ) < Pl \Q\p 2 >= i— < P1IQ5IP2 > 

m n 

= i{27if5{p 1 -p 2 ){Ei - E 2 -m n ) < Pl \ ^ 1 > 2 > . (48) 
If one has Ei 7^ E 2 in the symmetry limit, then Eq. (jl8|) amounts to : 

< Pi\Q\P2 >= z(2tt) 3 5(Pi - p 2 ) < Pl \ ^ K > 2 >, (49) 
which vanishes by the PCAC hypothesis. Therefore, if one evaluates the matrix element of 



[Q, A] between nucleon states according to 



< N(p')\[Q,A}\N(p) >= £ < W)\Q\X >< X\A\N(p) > -c.t, (50) 

x 

where the sum includes integration over the momenta and c.t. denotes the crossed term, 
then the only states which survive in the chiral limit are those for which : 

E x = E N (p) or E x = E N {p'). (51) 

This limits the sum over X to the following possibilities : 

1. X is the nucleon (Fig. Et)- 

2. X is the semi-disconnected nN state (Fig. Eb)with the soft pion anihilated by Q (the 
crossed term vanishes due to Eq. lEfEl) ). This is, up to a factor, the soft pion production 
amplitude we are looking for. 

3. X is the semi-disconnected 7tN state (Fig. Efc) with the soft pion created or anihilated 
by A 

Other intermediate states, for instance an isobar replacing the nucleon in Fig. Et or a heavy 
meson replacing the pion in Fig. Et, cannot be degenerate with the initial or final state 
and therefore are suppressed by PCAC. This is also the case for pions loops some of which 
are shown on Fig. El Here the ener gy g ap is due to the fact that the intermediate pions 
carry momentum. As is well known [23], the vanishing chiral limit of these contributions 
is in general reached in a non analytical way due to the infra-red part of the momentum 
integration. The leading non analytic term is due to the loop shown in Fig. Et and its 
infra-red part involves again the production of a soft pion by the operator A. So it could 
be computed within our framework without extra hypothesis. However, at this stage of 
investigation of DVCS, we think it is reasonnable to keep only the terms which do not 
vanish in the chiral limit, i.e. those shown on Fig. El 

We note that the term corresponding to Fig. Eh does not contribute to pion electro- 
production since the electromagnetic current cannot create a pion out of the vacuum. By 
contrast it can contribute to pion twist-2 production if the momentum of the intermediate 



3,® 
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Figure 5: Soft pion expansion of the commutator. Note that the intermediate lines are not 
propagators. 



p p p 

(a) (b) 
Figure 6: Example of pion loops. 

pion, that is q— <f in this case, is of order m„.. Since £ 7 = (q — q') 2 < 0, we have (q~q f ) 2 > |t 7 | 
and therefore we can neglect this term if we impose that |t 7 | be much larger than m^. More 
precisely we ask that £ 7 or i be a typical hadronic scale which need not be large but which 
does not vanishes when we let m n go to zero. Clearly this is compatible with the conditions 
(120)1 for factorisation of the DVCS reaction. Though this is not immediately necessary we 
from now on assume that this condition is realised independently of the reaction considered, 
that is whatever the operator A. 

Under this condition we can write the soft pion expansion of the commutator as : 

<N(p')\[Q a ,A}\N(p) >= 
dK 



(2tt)*2E n (K) 



(< N(p')\Q a \N(K) >< N{K)\A\N{p) > -at.) + 



/ 



dK 



(2n) 3 2E 7T (K) 



< 0\Q a \ir p (K) >< N(p')ir IJ (K)\A\N(p) > +0(m*). (52) 



where the sum over the spin and isospin labels is understood. From the definition of the 
isovector axial form factors : 



< AV)|g 7 M 75yn u( M0)|iV(p) >= u(p' 



9A{t)-f + h A (t) 



ip' — pY 

2M 



7fef«(p) (53) 



we have : 



< iV(p')|Q Q |iV(p) >= (27r) 3 5(p / -p)^(0) M t (p) 7 5^(p) 



which, together with Eq. (|T7|l . leads to the following expression of the commutator 



<N(p')\[Q«,A]\N(p)> = g A (0) 

- 9a(0) 



^TyJ^Vhsf u(p') < N(p')\A\N(p) > 
1 < N(p')\A\N(p) > tt t (p)75y«(p) 



2E N (p) 

+ if* < N(p')n a (m n ,d)\A\N(p) > +0(m n ] 
Let us now restore the spin label and define the on mass shell vertex T by : 



(54) 



(55) 



< N(K',a')\A\N(K,a) >= u{K' ,a')T{[K'},[K})u{K,a) (56) 
where [K] is our special notation for the on shell 4-momentum that is 



[K] = K, [K} = \jK 2 + M 2 . 



(57) 



It is then a trivial task to show that Eq. (|55|) is, up to terms of order m w , equivalent to : 



< N(p')\[Q a , A}\N(p) >= -iUT B + if„ < N(p')7r a (m n ,0)\A\N(p) > +0(m n ), (58) 



with the Born term T B defined as : 



2U 



^.(p' + k^-M 
1 



T(\p' + K],\p\) 



r(p-K)-M 



u(p), 



(59) 



and k n = (m n , 0) is the 4-momentum of the pion at rest. The expression is formally 
covariant so we can use it in a frame where the pion is not at rest. The 0(m n ) corrections 
to the soft pion expansion then become 0(k n ) where k n is the 4-momentum of the moving 
pion. In this way one generates correctly only the terms which go like the velocity (~ k^/m^ 
) of the pion. The other terms of order A;^ in (|59|) cannot be distinguished from the other 
corrections to the soft pion expansion. 

Note that the expression of the Born term Tb is not exactly what one would expect from 
an effective Feynman diagram with pseudo-vector nN coupling because, in this case, the 
vertex Y would appear in the form T(p' + k n ,p) or T{p',p — k n ). This mismatch is due to 
the fact that in the soft pion expansion the vertex F comes naturally into play with its 
arguments on the mass shell, and therefore cannot conserve energy. In practice this energy 
non conservation is of order k n so one could replace [p' + k n ] and \p — k n ] by their off mass 
shell values p' + k n and p — k n since this would amount to change the 0(k n ) correction. 
However the advantage would be only of cosmetic nature because we do not know what is 
T when its arguments are off the mass shell! So we shall retain expression (|59|) . 

B. Evaluation of the commutators 

In the symetry limit the axial charge is time independent and we now that the symetry 
breaking term Q 5 vanishes as m% . So to compute the matrix element of the commutator : 

< N\[Q,A]\N >=< N\[Q 5 + —Q 5 ,A]\N >, (60) 

m n 

one can use the commutations rules : 

[Q%,q(x)] = - 7s ^n vd q(x), (61) 

and neglect the term Q 5 /m n . The error will be of order m n unless there is a kinematical 
enhancement due to the proximity of a pion pole in the matrix element (j6fljl . To check this 
we look for such poles in the amplitudes and the dangerous ones are shown on Fig. In 
the case of electroproduction the pion pole can only be in the NN channel as shown on Fig. 
UK- In the case of associated DVCS it can appear in the 77 channel (Fig. 0b), in the NN 



(a) 



(b) 





(c) (d) 
Figure 7: Pion poles in the amplitudes for electro-production (a) and associated DVCS (b, c, d). 

channel (Fig. Efc) or in the (771) channel (Fig. 01). The 77? channel is suppressed by powers 
of 1/Q , so we are left with pion poles which are of the form : 

1 1 

or 



t — m\ t 7 — ml 



which are not dangerous since we have assumed that t or t 7 remain finite in the chiral 
limit. 

Using Eq. (|61|) it is straightforward to compute the relevant commutators. From the 
definitions (fTT| l33|) one gets : 



f 



■ u/3 



(62) 
(63) 

(64) 



where Sy a are the iso- vector twist-2 operators introduced in (|32l33j) . 



C. Soft pion theorem for the pion electroproduction 



From Eqs. (|58|62|) we get the soft pion theorem, which we write directly for a moving 
pion : 

< NW^iKWlNb) >= T B + ±re aW < N(p')\jf\N(p) > +0(K), (65) 

Jit 



where the Born term is obtained from Eq. ()59|) with : 

r([K'],[K}) = [Ff([t])A^ + mt])^^)Y 



+ i (mm 1 -^ + ™^ ^ ([J ?,7 W) » , (66) 



and [t] = {[K 1 
Nambu et al. 



- \K 

21M 



,2 



2M 



This is nothing but the low energy theorem originally derived by 



. Using our assumption that t is non zero in the soft pion limit, one 



can check that the amplitude (|63| respects gauge invariance in the form 

(p' + K- P)„ (t b + je aW < N(p')\jf\N(p) > +0(kSj = 0, (67) 

that is the additional terms needed to have exact gauge invariance are in the 0(k n ) correc- 
tions The axial current matrix element which appears in Eq. (|65|) is given by Eq. (|53f) 
and involves the form factors su(^) and ^a(*)- The latter is not so well known experimentally 
but this does not matter because it is multiplied by p' — p = — k w where k^ = q — q' is 
the photon exchanged in the BH process. In the Lorentz gauge the term proportional to k^ 
does not contribute to the amplitude and, because t does not vanish in the soft pion limit, 
the term linear in k n can be pushed in the 0(k n ) corrections despite the pion pole in h,A(t). 
The rest of the calculation is just algebraic manipulation of expressions (|65I66|) so we do not 
need to give the details. For completeness we just mention the relation between the charged 
and cartesian pion production amplitudes : 

T(n±) = nnl ^ 2 \ T(^)=T(^). (68) 



D. Soft pion theorem for the pion twist-2 production 



The soft pion theorems for pion twist-2 production folllow from Eqs (|58|63l64|) : 

< N(p')n a (K)\J2Q 2 f S f (x,n)\N(p) >= 

f 

T B + -U^ < N(p')\S 5 vp (x,n)\N(p) > +0(k n ), (69) 

< N{p')ix a {K)\Y^Q)S%x,n)\N{p) >= 

f 

T% + -U*^ < N(p')\S V p(x,n)\N(p) > +0(k n ), (70) 

& J TV 

where the Born terms T B ^ are obtained from Eq. ([59|) with the vertices T^ 5 ) defined by (see 
Eq.®) : 

<N(K')\J2QfS f (x,n)\N(K) > = u(K')T([K'},[K},x,n)u(K), (71) 
/ 

<N(K')\J2Q}S%x,n)\N(K) > = u(K')T 5 ([Kl[K],x,n)u(K). (72) 
/ 

Note that we have restored the explicit dependance on x, h. Using the parametrisation (|36|) 
we have : 

r([K'],[K],x,n) 



([K] + [K']).n 



E Q} (h,M[x], [£], [A] 2 h.n + iE f/N ([x], [£], [A] 2 ) a ^ A]l/ ) , (73) 
T 5 ([K'],[K},x,n 



([K] + [K']).n 

E Q} (h,/n([x], [£], [A] 2 ) 7 .n + E f/N ([x], [£], [A] 2 )^) 7 5 , (74) 



where the bracketed variables are, by definition (see Section |I 



From Eq.dSnj we see that one has either ([K] = \p), [K'\ = [p' + K]) or ([K] = \p~K], [K'\ = 
\p'))- It is then apparent that the approximation ([x] —> x [£] — > £, [A] —* A) in Eqs (|73|74l) 



(5) 

only induces corrections of order k n . Similarly the matrix elements of Sy which appear in 
Eqs. (j69|70|) can be written as : 



< N{p')\S Va {x,n)\N(p) >-- 



(p + p').n 



x 



u(p') 



H v (x, i, t)^.n + iE v (x, i, t) 



a^n^p' - p) t 
2M 



T N 

-fu{p), 



(75) 



< N(p')\S Va (K,h)\N(p) >-- 



x 



u(p') 



Hv(x, f , t)j.n + E v (x, f , t) 



(p + p').n 
n. (p' — p) 



2M 



-fu(j>), 



(76) 



where now the hatted variables are : 



x 



X 



(p + p').n' 



(p' — p).n 
(p + p').h' 



and t = (p' — p) 2 . Again the approximation x —>■ x, £ — > £, (p' — p) — > A induces only 
corrections of order k n . This completes our derivation of the soft pion theorem for twist-2 
production. 

A few comments are in order. First we see that, as promised in the introduction, the soft 
pion theorem for twist-2 production involves no new quantities with respect to the elastic 
case. Second we note that the commutator term, which involves only the isovector operators 

(5) 

S v , depends essentially on the valence quarks while the Born term, which sums over all 
favours, is sensitive both to the valence and the sea quarks. Third we recall that gauge 
invariance is not an issue for the pion twist-2 production. The couplings to the photon fields 
have been factorized from the beginning, see Eq.JSHJ), independently of the final hadronic 
state. As in the elastic case, gauge invariance can be restored by a term which is explicitely 
of higher order in t/Q , and thus does not contribute in the Bjorken limit ["b^ . 



VI. ESTIMATE OF THE A(1232) CONTRIBUTION 

The soft pion theorem proposed in the previous section can be valid only for small values 
of the pion momentum. One can expect strong deviations when the mass of the pion-nucleon 
system approaches the first nucleon excitation, that is the A(1232). So to improve the range 
of validity of our estimates we now calculate the contribution of this resonance. As a starting 



N(P) A(P A ) 
Figure 8: Delta excitation by the twist- 2 operator. 



point, we consider the matrix element as shown in Fig. |Hl where the A-resonance is on its 
mass-shell. Subsequently, we discuss the modification due to the A — ► nN strong decay. 
The N — > A matrix element for the vector twist-2 operator has the following parametriza- 
3,Q Q : 



tion 



<A(p A )\J2Q 2 f S f (x, n)\N(p)> 
f 

\¥{p±)Tt {H M (x,Z,A 2 ) (-fcjj) + H E (x,£,A 2 ) {-K 



+ H c (x,Z,A 2 ) (-Kg,) ^} u(p), (77) 

where ^(pa) is the Rarita-Schwinger spinor for the A-field and T\ is the isospin 1/2 — > 3/2 
transition operator, satisfying : 

(\t* I T{ | \t z ) = (If, , 1A | \t.) , (78) 

with t z (T z ) the isospin projections of N (A) respectively, and where A = 0, ±1 indicate 
the spherical components of the isovector operator . The factor 1/6 in Eq. ((77jl results 
from the quadratic quark charge combination (e^ — e 2 d )/2. Furthermore, in Eq. the 
covariants K,%] E ' C are the magnetic dipole, electric quadrupole, and Coulomb quadrupole 



covariants 



r u 3(M A + M) A 

2M((M A + M) 2 - A 2 ) ^ Act 

^ - MZ(A 2 ) £ ^ A <° 7 ' 1 



where we introduced the notation : 

Z(A 2 ) = [(M A + Mf - A 2 ] [(M A - M) 2 - A 2 ]. (80) 

Here P = (pa + £>)/2, A = pa _ P 5 and p A = M A (Ma = 1232 MeV). As we are considering 
the production of a A here, we explicitely adopt the notation A 2 in all expressions in this 
section to denote the momentum transfer squared to the hadronic system, in order not to 
confuse with the variable t introduced before. In Eq. ifTTjl. the GPDs H M , H E , and Hq for 
the N — > A vector transition are linked with the three N — > A vector current transition 
form factors G* M , G* E , and G* c through the sum rules : 

dx H M ^c(x, £, A 2 ) = 2 G* M>E>C {A 2 ) , (81) 

where G* MEC (A 2 ) are the standard magnetic dipole, electric quadrupole and Couloub 
quadrupole transition form factors respectively [26]. As is well known from experiment, 
the N —>■ A vector transition at small and intermediate momentum transfers is largely dom- 
inated by the N — ► A magnetic dipole excitation, parametrized by G* M (A 2 ). At A 2 = 0, 
its value extracted from pion photoproduction experiments is given by G* M (0) ~ 3.02 j^f. 
For its A 2 -dependence, we use a recent phenomenological parametrization |2j| from a fit to 
pion electroproduction data. Given the smallness of the electric and Coulomb quadrupole 
iV — > A transitions, we will neglect in the following the contribution of the GPDs H E and 
H c . 

The N — > A matrix element for the axial twist-2 operator has the following parametriza- 
tion 3,Q : 



<A(p A )\J2Qpf(x, n)\N(p)> 
f 

1 ^(pa) Tt sf- {Ci(x,£, A 2 )h[s + C 2 (x,£, A 2 ) A? A ' " 



+ C 3 (x, C, A 2 ) — (7. A np - 7 .n A p ) 



+ C,(x, £, A 2 ) A (A . - A^)} u(p) . (82) 



The GPDs Ci,C 2 ,C 3 and C 4 entering in Eq. (|82j) are linked with the four A" — > A axial- 
vector current transition form factors C^(A 2 ), Cq(A 2 ), C^(A 2 ) and C^(A 2 ), introduced by 
Adler (3| through the sum rules : 



j dxCi 


>,e,A 2 ) 


= 2 C 5 A (A 2 


j ^dxC 2 


>,e,A 2 ) 


= 2 C 6 A (A 2 


j dxC% 


>,e,A 2 ) 


= 2 C 3 A (A 2 


j dxC\ 


>,e,A 2 ) 


= 2 C^A 2 



(83) 

For small momentum transfers A 2 , PCAC leads to a dominance of the form factors C§ and 
Cq. For C§, a Goldberger-Treiman relation for the A" — > A transition leads to : 

Jlc*W = -£f L 9A- (84) 

Using the phenomenological values f n NA — 1-95, f W NN — 1-00, and qa — 1.267 one obtains 
C^(0) ~ 1.01. The form factor C^(A 2 ) at small values of A 2 is dominated by the pion-pole 
contribution, given by : 

M 2 

Of (A 2 ) = ——-cm- (85) 

n 

Because we are only interested in the limit —t, — 1 7 >> m 2 , we neglect the pion pole contri- 
bution of C 6 in the following, consistent with the discussion of Section IV Bl 

For the two remaining Adler form factors C% and Cf, a detailed comparison with ex- 
perimental data for //-induced A ++ production led to the values [3] : C^(0) ~ 0.0 and 
C^(0) ~ —0.3. Given the smallness of these values, we will neglect in the following the 
contributions from the GPDs C3 and C4. 

To provide numerical estimates for the A" — > A DVCS amplitudes, we need a model for 
the three remaining 'large' GPDs which appear in Eqs. 11771821) i .e. H M , C\ and C 2 . Here 
we will be guided by the large A" c relations discussed in lj. Iis|. 



These relations connect 



the iV -> A GPDs H M , d, and C 2 to the N ^ N isovector GPDs E u — E d , H u — H d , and 
E u — E d respectively. These relations are given by 0] : 



Ci(x,£,A 2 ) 
C 2 (x,£,A 2 ) 



2 

73 
v/3 

V3 
4 



E u (x,^, A 2 ) 

^(x,e,A 2 ) 
^(x,e,A 2 ) 



^(x,e,A 2 ; 

iJ d (x,e,A 2 ; 

^(x,e,A 2 ; 



(86) 



To give an idea of the accuracy of these relations, we calculate the first moment of both 
sides of Eq. (|86f and compares their values at A 2 = 0. For H M . we obtain for the Ihs 
the phenomenological value 2G* M (0) ~ 6.04, in comparison with the large N c prediction of 
(2/V3)k v ~ 4.27, accurate at the 30 % level Q- 

For Ci, the phenomenological value yields 2C 5 (0) ~ 2.02, whereas the large N c prediction 
yields v^SU — 2.19, accurate at the 10% level. For the pion-pole contribution to C 6 , we 
obtain the same accuracy as for C 5 . Furthermore, for the N —> A DVCS in the near forward 
direction, unlike the N — > N DVCS case, the axial transition (proportional to the GPDs 
C\ ) is numerically more important than the vector transition (parametrized by Hm). This 
is because Hm is accompanied by a momentum transfer A in the tensor K M as is seen 
from Eq. (f79|) . in contrast with the structure associated with the GPD C\. From these 
considerations, we estimate that the large N c relations of Eq. (|86|) for the N — > A transition 
allow us to estimate at the j*P ~^ 7^ DVCS process, at the 30 % accuracy level or better. 

In order to calculate the coherent sum of the non-resonant ADVCS process discussed 
before and the A-resonant process, we next discuss the modification of the 7*iV — ► 7A 
matrix elements due to the A — > nN strong decay. 

The matrix element of the vector operator for the *y*N — > 7 A transition followed by 
A — > nN decay, is modified from Eq. (JT7J) to : 

<7iN\J2Q}S f (x, n)\N> Delta 
f 

-fnNA i(j-PA + W) 



-1- 



x 



> W 2 _ Ml + iwr A (w) 

1 , , s , , , 2 



9aP ~ g7a7/3 



3W 



MPa)/3 - lp{PA)a) ~ 7^(Pa)«(Pa)i3 



-K 



M \ @ K ~ 



u(p), 



(87) 



where p A = p' + k w , A = p' + k n — p, and where we only indicated the leading transition, 
proportional to Hm- The isospin factor X takes on the values : X = 2/3 for the n°p final 
state, and X = —^/2/3 for the n + n final state. In the A propagator in Eq. (18*71) . the energy- 
dependent width T A (W) is given by : 



r A {w) = r A (M A ) 



K{w) 



Ma 
W ' 



(88) 



with T A (M A ) ~ 0.120 GeV, and where \k* n (W) \ is given in Eq. ijTTH . 

Analogously, the matrix element of the axial- vector operator for the 7*iV — > 7A transition 
followed by A — > rrN decay, is modified from Eq. (jH2j) to : 



<7riV|^Qj5j(a;, n)\N > Delta 
f 

-Una a — /„a * (7-Pa + W) 



X- 



K a u(p 



W 2 - Ml + iWT A (W) 



x 



x 



11 2 

£a/3 - g7a7/3 - 7^ (7a(PA)^ - lp{.VA) a ) ~ 7^{.PA) a {.PA)l3 



(89) 



where we again only indicated the leading transition proportional to the GPD C\. 

Finally, in order to calculate the YP ~^ •yn'N process, we also need to specify the BH 
process of Eq. (j2J) associated with the N — > A transition. The corresponding A contribution 
to the matrix element J v entering in Eq. (J2J) is given by : 



<nN\f(0)\N> Delta 



W 2 - Ml + iWT A (W) 



g aP - 37.7/3 - ^ 



{ja(j>^)p - 7/3 (PA) a) 



3W 2 



(pa) q (pa)/3 



x{ G* A/ (A 2 ) (-K"f + G* E (A 2 ) {-K E f + Q(A 2 ) (-/C c f} u(p), 



(90) 



in terms of the same vector N — > A transition form factors G* M E c as discussed before. The 



isospin factor X is the same as specified following Eq. lj%TJ) . 



VII. RESULTS 

In the absence of available data for the 7*p — > jtiN process, we can get an idea on the 
accuracy of our estimates by comparing the pion production amplitudes of Eqs. (|65|90j) . 
which enter in the Bethe-Heitler process, with the pion photo- and electroproduction cross 
sections. We are interested here in the region of not too large — £ 7 < 1 GeV 2 , corresponding 
with the virtuality of the photon in the Bethe-Heitler process. In this kinematical range, a 
large amount of pion photo- and electroproduction data exist to compare with. In particular, 
the pion photoproduction cross section ■jp — > nN (proportional to the cross section of the 
associated Bethe-Heitler (ABH) process for £ 7 — > 0), is given by 

(£)„ = W^f ^ \ E I «.(., ^ < ^m\N > |». (9!) 

where is the photon cm. energy, e u (q, A) is the photon polarization vector, W^at is 

the 7riV cm. energy, and f is the electromagnetic current operator. We show the results 
of the different model contributions discussed above to the pion photoproduction total cross 
sections in Fig. EH As can be noticed from Fig. El our estimates consisting of a soft-pion 
production amplitude supplemented by a A-resonance production mechanism reproduce the 
cross sections on the lower energy side of the A(1232) resonance. Around resonance po- 
sition our simple model overestimates the cross sections by about 10 %, which is mainly 
due to rescattering contributions which can be included by a proper unitarization of the 
amplitude, which we did not perform in our simple estimate. At the higher energy side of 
the A(1232) resonance, the present model somewhat overestimates the data. Besides the 
unitarization, this due to the increasingly important role from the t-channel exchanges of p 
and ui vector mesons in the non-resonant part of the amplitude. It is known that these vector 
meson exchanges yield a destructive interference on the higher energy side of the A(1232) 
resonance. Since our objective here is not to present a phenomenological model for pion 
photoproduction, able to precisely describe the available data, but to provide an estimate 
for the associated DVCS process, we do not include a proper unitarization or vector meson 
exchanges as one is not able at this point to model the corresponding contributions for the 



ADVCS process. Therefore the quality of the description shown in Fig. El is indicative of the 
quality of our corresponding estimates for the ADVCS process which are performed along 
the same lines, i.e. by the sum of a non-resonant soft-pion production amplitude and a 
A(1232) resonance production amplitude as outlined above. 

In the following calculations we use, for the nucleon GPDs, the model described in 
Refs. 0,0,13, 

to which we refer the reader for details. We construct the GPDs from a 
double distribution based on the forward unpolarized quark distributions of MRST01 
(for the GPD H) and based on the forward polarized quark distributions of Ref. {^J (for the 
GPD H). To construct the double distribution, we use a profile function with parameters 
Kai = b sea = 1 as detailed in Ref. 3|. For the t-dependence of the GPD (at the moderately 
small values of t considered in this paper) we adopt, unless stated otherwise, a factorized 
ansatz by multiplying the x and £ dependent function by the corresponding form factor (in 
t) so as to satisfy the first sum rule. In the calculations for the beam charge asymmetry, we 
also compare the results with a model for the GPDs where a D-term is added to the double 
distribution, with the parametrization given in Ref. p3|. For all other calculations, where 
it is not stated explicitely, the results do not include a D-term in the model for the GPDs. 

In Fig. EH we study the different ADVCS processes and show their contributions to the 
7-fold differential e~p — > e~j7rN cross sections, differential in Q 2 , xg, £ 7 , <3>, the ttN invariant 
mass W n N, and the pion solid angle f2* in the ttN rest frame. By comparing Figs. [HI and 
ITflj one sees that the ratio of the non-resonant to resonant contributions is larger for the 
ADVCS process compared to the pion photoproduction process. 

In Fig. EJ we compare the 5-fold differential e~p — ► e^jirN cross sections, i.e. integrated 
over the pion solid angle fi* , for the ABH, ADVCS and ABH + ADVCS processes Q 
for JLab kinematics. Clearly, the ABH largely dominates the cross sections. The resulting 
beam spin asymmetries (BSA), for a polarized lepton beam, are around 5 - 10 % for the 
e~p — > e~jTT + n process. For the e~p — > e~77r°p process, the BSA grows when approaching 
the iiN threshold, where it reaches the same value as for the e~p — > e~jp process. This 
can be easily understood because at threshold, the amplitude for the ep — > ej^p process is 
obtained from the ep — > e^p process by attaching a soft pion to the initial and final proton. 
This is what we called the Born term in Section 5 (Fig. El a). This amounts to multiply the 
DVCS and BH amplitudes of the ep — > ejp process by the same factor when calculating 
their ep — > e^n°p counterparts. Therefore, when taking the ratio of cross sections in the 
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Figure 9: Total pion photo-production cross sections for the different model contributions considered 
in this paper. Dashed curve : commutator contribution. Dashed-dotted curves : commutator + 
Born contributions. Dotted curves : A contribution. Solid curves : commutator + Born + A 
contributions. The data are from Ref. (diamonds), Ref. |ilj (circles), and Ref. ji^] (triangles). 

BSA, which is due to the interference of ABH and ADVCS, this common factor drops out 
and one obtains the same BSA as for the e~p — > e~7p process. Note that this is not the 
case for the e~p — > e~'j7r + n process, where both commutator and Born terms contribute. 
Futhermore in the Born term for the p — > ir + n transition, amplitudes involving both proton 
and neutron GPDs interfere according to whether the charged pion is emitted from the final 
or intial nucleons respectively. This results in a much smaller BSA at threshold for charged 
as compared to neutral pion production. When moving to higher values of W n N the A 
contribution becomes important, and the ratio of ADVCS to ABH changes compared to the 
e~p — > e~7p process. Around A(1232) resonance position, the BSA for the charged and 
neutral pion production channels reach comparable values, around 10 %. 

If one does not perform a fully exclusive DVCS experiment, one actually measures the 
cross section 



da exp 



da{ep — > e^p) 
dQ 2 dx b dt <i$ dQ 2 dx b dt d$ 



(1 + R 



inel i 



(92) 




Figure 10: Different contributions to the 7-fold differential cross sections for the ADVCS processes in 
JLab kinematics, with pion emitted in the same direction as the recoiling nN system (corresponding 
with 9* = 0°). Dashed-curve : commutator contribution. Dashed-dotted curves : commutator + 
Born contributions. Dotted curves : A contribution. Solid curves : commutator + Born + A 
contributions. For the N — > A GPDs both C% and Hm are included. 

with Ri n ei the ratio of the integrated inelastic ep — > cyirN cross section to the cross section 
for the ep — * ejp reaction (i.e. the 'elastic' DVCS process), given by 

/ f w ™* dw da(ep -> e-firN) \ da(ep -> e^p) 
\JM+m„ nN dQ 2 dx b dt (2$ dW nN J dQ 2 dx B dtd§ ' 

This is the ratio k introduced in Eq. (|19|) . but written in a more familiar way. The ratio 
Rimi depends upon the upper integration limit W max , determined by the resolution of the 
experiment. 
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Figure 11: 5-fold differential cross sections (upper panels) and corresponding beam-spin asymme- 
tries (lower panels) for the e~p — > e~jirN reactions in JLab kinematics. Dotted curves : ABH; 
dashed curves : ADVCS; solid curves : ABH + ADVCS. The arrow gives the elastic value of the 
BSA for the BH + DVCS process, corresponding with W wN = M = 0.939 GeV. 



We can now provide an estimate of the 'contamination' to the BSA for a not fully exclusive 
experiment. In an experiment where one does not separate the 7p and 771 N final states, one 
actually measures 



(BSA) exp = 



Ada el (1 + AR 



inel J 



2 da* 1 (1 + R inel ) ' 



(94) 



where Ada el and da el stand for 

Adad = dQ Z h xsdu^ ep -> eip) - d Q ZT B du^ ep -> e7p) ' 

^ = 5 { ^sdtd^ ep e7p) + WMd^ {ep e7p) ) ' (95) 

with /i the lepton beam helicity. In Eq. (jMj) . -Rj ne ; is given as in Eq. (|93|) . and AR inel is the 
corresponding ratio of inelastic to elastic DVCS helicity cross sections : 

AjR . = ( f Wmax dw Adajep -» e77rAQ \ Adajep -» e7p) 
\/M+m* 7tN dQ 2 dx B dtd^dW 7TN J dQ 2 dx B dtd§ 

From Eq. (JM|), one sees that for a not fully exclusive experiment the 'elastic' beam-spin 
asymmetry (BSA) el for the ep — > e7p process is related to the measured beam-spin asym- 
metry (BSA) exp through : 

Ada el 

(BSA) el = = Rbsa • (i?5A) e ^, (97) 

where the correction factor -Rbsa is given by : 

_Rg Sj4 = - — ^ inel (gg) 

1 + ARi ne i 



In Fig. H21 we show the ratios Ri ne i and AR ine i- One sees that when integrating up to 
W m ax — 1-35 GeV, the helicity cross sections ratio AR ine i reaches about 10 % for •yn+n 
and 77r°p final states separately. On the other hand, the unpolarized cross section ratio 
Rinei is much larger and reaches about 40 % for each channel separately. This difference 
originates in the different ratio of the ABH cross section as compared to the corresponding 
ADVCS ratio. This different ratio for the ABH and ADVCS processes has as consequence 
that the BSA, which is due to an interference of both, receives an important correction in 
an experiment where the final state cannot be fully resolved. 

In Fig. [TBI we show the resulting correction factor for the BSA defined in Eq. for JLab 
kinematics. For an experiment which measures an e~ and a proton, and which reconstructs 
the final 7 from the missing mass, but where the resolution does not permit to fully separate 
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Figure 12: Ratios R ine i and A.R inei of e"p -> e'jirN (ABH + ADVCS) to e"p -> e"7p (BH + 
DVCS) cross sections according to Eqs. as function of the upper integration limit W max in 

JLab kinematics. Dashed (solid) curves : ratio of unpolarized (polarized) cross sections respectively. 

the 7» final state from a 77r°» final state, such as in the bulk of the events of the first DVCS 

Q 

experiment at CLAS one obtains a correction factor of around 1.3 when integrating up 
to W max ~ 1.35 GeV in the kinematics considered in Fig. UHl For an experiment which 
only measures an e~ and a 7 in the final state and where the resolution does not permit to 
separate the p hadronic final state from the n + n and 7r°p hadronic final states, the correction 



factor, integrated up to W n 



1.35 GeV, amounts to about 1.6 in the same kinematics. 
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Figure 13: Correction factor according to Eq. to apply to the measured BSA to extract the 
e~p — > e~7p BSA, when the measurement is not fully exclusive. The correction factors are given 
for the situations where the ~/p final state is contaminated either by 77r + n (upper panel), by 77r°p 
(middle panel), or by both 77r + n and 77r°p (lower panel). The correction is given as function of the 
upper integration limit W max in JLab kinematics. 

In Figs. HU - HBl we study the corresponding effects in the kinematics accessible in the 
HERMES experiment. 

One sees from Fig. [TU that in the kinematics accessible at HERMES, the ABH still 
dominates the cross sections. The interference of the ABH and ADVCS processes leads 
to a BSA which is around 10% for the e~p — > e^jir + n reaction. For the e~p — > e~77r°p 
reaction, the BSA rises towards nN threshold where it reaches the value of the BSA for the 



e~p — > e _ 7p reaction, as discussed before. Around the A(1232) resonance position, the BSA 
for the e~p — > e~ r yiT p reaction reaches about 15 %. 

For the present DVCS experiments at HERMES [f| where the experimental resolution 
does not allow to fully reconstruct the final state, it is important to estimate the contribution 
of the 77r + n and •y7i°p final states, which we show in Fig. [13 When integrating the cross 
sections up to W max ~ 1.35 GeV, the helicity cross sections ratio AR ine i reaches about 3 % 
(5 %) for the •yiT + n (77r°p) final states respectively. On the other hand, the unpolarized cross 
section ratio Ri ne i reaches about 10 % for each channel separately. This different ratio leads 
to a correction for the BSA in an experiment where the final state cannot be fully resolved, 
which is shown in Fie;. For kinematics close to the first DVCS experiment at HERMES 

Q 

|o] which measured an e and a 7 and where the resolution did not permit to separate the p 
hadronic final state from the n + n and 7i°p hadronic final states, one sees that the correction 
factor on the BSA due to the A resonance region, i.e. integrated up to W max ~ 1.35 GeV, 
amounts to about 1.1 . 

We next discuss the beam charge asymmetries (BCA) between the e + p — > e +r y7iN and 
e~p — * e~jTrN processes. The BCA in kinematics accessible at HERMES is shown in Fig.lTTI 
for two models of the GPDs, one including the D-term and one without the D-term contri- 
bution. As for the BSA, one sees that for the j7i°p final state, the BCA reaches the same 
value as for the elastic DVCS process when approaching the nN threshold. Furthermore, 
one sees that since the D-term only contributes to the Born terms, it mainly manifests itsef 
in the neutral pion production channel around threshold, while its effect is very small on 
the charged pion production channel. Around A(1232) resonance, the effect of the D-term 
is negligible. 

The BCA between the e + p — > e +r yX and e~p — > e~^X processes has been measured at 
HERMES jsfj. Because this experiment does not allow to distinguish the hadronic final 
states X = p from X = n + n and X = n°p, it is of interest to estimate the 'contamination' 
by the associated pion production. An experiment which does not separates X = p from 
X = ttN measures : 

{BCA) ~ af + (l + R + )+a- (1 + R_y ( " } 
where o± stands for the cross section da(e ± p — > e ± ^p) / {dQ 2 dxs dt d$) of the 'elastic' 
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Figure 14: 5-fold differential cross sections (upper panels) and corresponding beam-spin asymme- 
tries (lower panels) for the e~p — > e~^TrN reactions in HERMES kinematics. Dotted curves : ABH; 
dashed curves : ADVCS; solid curves : ABH + ADVCS. The arrow gives the elastic value of the 
BSA for the BH + DVCS process, corresponding with W nN = M = 0.939 GeV. 

process, and where the ratios R± stand for : 



rWmax da(e ± p — > e ± jnN) 

(?f± Jm+ttit, w dQ 2 dxs dt G?$ dW-xN 



(100) 



From Eq. one sees that for a not fully exclusive experiment the 'elastic' (BCA) el is 

obtained from the measured (BCA) exp through : 



(BCA) 



el 



a 



-el 

e+ 
el 



R 



BCA 



(BCA) 



exp 



(101) 
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Figure 15: Ratio of e~p -> e~jirN (ABH + ADVCS) to e"p -> e"7p (BH + DVCS) cross sections 
according to Eqs. (JHSHHEl) as function of the upper integration limit W max in HERMES kinematics. 
Dashed (solid) curves : ratio of unpolarized (polarized) cross sections respectively. 



where the correction factor Rbca is given by : 



Rbca 



1 + {af + R + + af-R-)/(a e e i + a?-) 
1 + (af + R + - af.R.)/{af + - af.\ 



(102) 



The correction factor Rbca is shown in Fig. EH for HERMES kinematics. For an exper- 
iment which does not separate the p hadronic final state from the n + n and n°p hadronic 
final states, the correction factor on the BCA, integrated up to W max ~ 1.35 GeV, amounts 
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Figure 16: Correction factor according to Eq. to apply to the measured BSA to extract the 
e~p — > e~7p BSA, when the measurement is not fully exclusive. The correction factors are given 
for the situations where the ~/p final state is contaminated either by 77r + n (upper panel), by 77r°p 
(middle panel), or by both 77r + n and 77r°p (lower panel). The correction is given as function of the 
upper integration limit W max for HERMES kinematics. 



to about 1.8 for the model without D-term and reaches around 1.1 for the model with D- 
term. The much smaller correction in the presence of the D-term can be understood as in 
this case the elastic BCA has the same sign and similar magnitude as the inelastic BCA 
around A(1232) resonance. On the other hand, in the absence of the D-term contribution, 
the elastic BCA is small and negative, yielding a significant different result from the positive 
BCA around A(1232) resonance. Therefore, the resulting correction factor is much larger 



for the GPD model without D-term. As the correction of the BCA can be sizeable according 
to the model for the GPDs, this clearly calls for a fully exclusive measurement to separate 
the different final states, in order to reliably extract information on the GPDs. Such an 
experiment is planned in the near future at HERMES using a recoil detector |3tT |. 




Figure 17: Beam charge asymmetries (BCA) for the ABH + ADVCS processes in HERMES kine- 
matics. The result is shown for 2 different models of the GPDs : without D-term (dashed curves) 
and with D-term (solid curves). The corresponding arrows give the elastic value of the BCA for 
the BH + DVCS process, corresponding with W^n = M = 0.939 GeV. 
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Figure 18: Correction factor according to Eq. (|1()2|) to apply to the measured BCA to extract the 
BCA between the 'elastic' processes e^p — > e ± 7p, when the measurement is not fully exclusive. 
The correction factors are given for the situations where the ^p final state is contaminated either 
by 77r + n (upper panel), by 77r°p (middle panel), or by both 77r + ra and 77r°p (lower panel). The 
correction is given as function of the upper integration limit W max for HERMES kinematics. The 
solid and dashed curves correspond with the two GPD models as described in Fig. El 

In Figs. HHH221 we show the results for kinematics accessible at COMPASS Fig. ITTH 
shows the differential cross section and BSA for the [i + p — > fi + 'yrrN processes. In contrast 
to the previous results for JLab and HERMES kinematics, we see that at COMPASS the 
H + p — > fi+jiiN cross section is dominated by the ADVCS process. The interference with the 
small BH yields only a small value for the BSA. Due to the dominance of the ADVCS over 



the BH, one sees in Fig. [201 an opposite trend for the ratios Ri ne i and AR ine i of Eqs. (j93l l96|) 
in comparison with the ones shown in Figs. [12] and HH1 for JLab and HERMES respectively. 
The larger value of AR ine i compared to Ri ne i, in particular for p + p — ► p + jii p, leads to a 
correction factor Rbsa of Eq. (|HK1) which is slightly smaller than one. For an experiment 
which does not permit to separate the p hadronic final state from the n + n and n°p hadronic 
final states, the correction factor on the BSA due to the A resonance region, i.e. integrated 
up to W max ~ 1.35 GeV, amounts to a value around 0.95 at COMPASS. 

In Figs. Ell E2 we study the BCA at COMPASS. To obtain sizeable interferences with 
the BH process, we show the results for a lower beam energy of 100 GeV as also accessible at 
COMPASS. We compare the results for two models of the GPDs. The first model consists 
of a factorized ansatz for the t-dependence (compared to the x- and £ dependences) of the 
GPDs as used in the previous calculations, and does not include the D-term. The second 
model includes the D-term and uses an unfactorized Regge ansatz for the GPDs as specified 
in Ref. As can be seen from Fig. EH the it production process is mainly sensitive to 

the differences between those models and at the threshold the BCA for the p ± p — > p ± 'yii p 
and p ± p — > p ± jp reactions are the same. For an experiment which does not separate the p 
hadronic final state from the 7r + n and 7r°p hadronic final states, we estimate the correction 
factors on the BCA in Fig.[22j When integrating the nN spectrum up to W max ~ 1.35 GeV, 
the correction factor on the BCA asymmetry due to both n°p and ir + n channels amounts 
to about 1.35 for the model without D-term and reaches around 1.05 for the model with 
D-term. 
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Figure 19: 5-fold differential cross sections (upper panels) and corresponding beam-spin asymme- 
tries (lower panels) for the /j, + p — ► [x + ^itN reactions in COMPASS kinematics. Dotted curves : 
ABH; dashed curves : ADVCS; solid curves : ABH + ADVCS. The arrow gives the elastic value of 
the BSA for the BH + DVCS process, corresponding with W nN = M = 0.939 GeV. 

VIII. CONCLUSION 



We have developed a model to calculate the cross section for producing an extra low 
energy pion in the photon electro-production reaction. Our primary goal is to provide a 
reasonable estimate of the contamination of the 'elastic' process by this associated reaction 
when the experimental data are not fully exclusive. For the various observables which 
are generally considered of interest we have defined correction factors by integrating the 
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Figure 20: Ratio of fi + p -> fi+^nN (ABH + ADVCS) to fi+p [i + jp (BH + DVCS) cross sections 
according to Eqs. (jHSHHEll as function of the upper integration limit W max in COMPASS kinematics. 
Dashed (solid) curves : ratio of unpolarized (polarized) cross sections respectively. 



associated reaction cross sections up to a given cutoff. 

To build our model we have used the time honored soft pion technique based on current 
algebra and chiral symmetry. In the case of the DVCS reaction, which we always consider 
in the Bjorken limit, we have assumed that it was possible to first invoke the factorization 
theorem and then to use chiral symmetry to evaluate the matrix elements of the twist 
two operators involving one soft pion. Our derivation applies to the kinematical region 
ml « — t « Q 2 of DVCS type processes, which corresponds with the kinematical range 




Figure 21: Beam charge asymmetries (BCA) for the ABH + ADVCS processes in COMPASS 
kinematics. The result is shown for 2 different models of the GPDs. Dashed curve : GPD model 
with factorized t-dependence and without D-term. Solid curve : GPD model with unfactorized 
t-dependence and with D-term. The corresponding arrows give the elastic value of the BCA for the 
BH + DVCS process, corresponding with W nN = M = 0.939 GeV. 

of experiments considered at JLab, HERMES and COMPASS. The order in which one 
applies the chiral limit (m n — > 0) and the Bjorken limit (Q 2 — ► oo) is a point which certainly 
deserves further attention, see e.g. Ref. |38l | where new low-energy theorems were derived 
for the 7*A^ — > ttN process at large virtualities of the photon (Q 2 » An^n/m^). In this 
respect, the chiral perturbation theory approach developed in Refs. (sol. l4fl| to study the 
quark mass dependence of the parton distribution may be useful. For such a systematic 
chiral perturbation expansion to converge, one also sees that the momentum transfer —t 
should be bound from above, and be smaller than the chiral symmetry breaking scale of 
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Figure 22: Correction factor according to Eq. (|1()2|) to apply to the measured BCA to extract the 
BCA between the 'elastic' processes fi^p — > /i ± 7P, when the measurement is not fully exclusive. 
The correction factors are given for the situations where the ■jp final state is contaminated either 
by 77r + n (upper panel), by 77r°p (middle panel), or by both jir + n and 77r°p (lower panel). The 
correction is given as function of the upper integration limit W max for COMPASS kinematics. The 
solid and dashed curves correspond with the two GPD models as described in Fig. I2"T1 

order (47r/ 7r ) 2 . 

The advantage of our approach is that the GPDs necessary to calculate the associated 
pion production are the same as in the elastic DVCS process. So, in a relative sense, our 
estimate is largely independent of the details of the GPDs. 

To extend our estimate to pions of higher energy we have added the P-wave production 



assuming it is dominated by the A(1232) isobar production. Here, following the approach 
of Ref. (l2j we have used the large N c limit to relate the GPDs of the N — ► A transition to 
the N —>■ N ones. Again this minimizes the model dependence of our results. 

For those experiments which do have the resolution to measure the YP l^N process, 
in particular in the A(1232) resonance region [ijj, our calculation provides a prediction using 
the large N c limit for the N — > A GPDs. The measurement of the YP ~~ * 7 A — ► l^N process 
holds the prospect to access information on the quark distributions in the A resonance, which 
are totally unknown at present. 

Generally the DVCS amplitude must be added coherently to the Bethe-Heitler amplitude. 
This is also true for the associated reactions and we have used the same approach to estimate 
the associated BH and DVCS amplitudes. For the BH amplitude this amounts to calculate 
the pion electro-production amplitude and this gives us a chance to rate the validity of our 
results by comparing to the existing data. We find that, for the integrated cross sections, 
our estimate are probably valid up to W max = 1.35 GeV. 

We have performed our calculations for a set of kinematical conditions which are repre- 
sentative of the present or planned experiments at JLab, HERMES and COMPASS. In a 
regime where the ADVCS process dominates the cross section, such as is the case at COM- 
PASS, we find that the pionic contamination (integrated up to W max = 1.35 GeV) never 
exceeds 10%. In a kinematical regime where the ABH process dominates, such as is typically 
the case at HERMES and in particular at JLab, the pionic contamination may become much 
larger and calls for fully exclusive experiments. In particular, it has an effect on the beam 
spin and beam charge asymmetries. For instance the correction to the BSA due to charged 
(neutral) pions in JLab kinematics can reach 30% each for an experiment which is not able 
to distinguish a N from a tcN for a cutoff W max = 1.35 GeV. The effect on the BSA due to 
the n°p and n + n production in HERMES kinematics is of the order of 10%. On the other 
hand the correction factor on the BCA to obtain the 'elastic' BCA from an experiment not 
able to distinguish a N from a nN can be as large as a factor 1.8 at HERMES and 1.35 at 
COMPASS depending on the model for the GPDs. 

To summarize, for a cutoff of W max = 1.35 GeV the correction to the cross sections is 
moderate but for the BSA and BCA our calculations indicate that it is wise to consider fully 
exclusive experiments. 
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Strictly speaking, exclusivity then means that there is one photon and one proton in the final 
state and nothing else. In practice the final state always contains soft photons which build the 
radiative tail but their effect is included in the calculable radiative corrections. So we ignore 
this subtility here and take Eq.Q stricto sensu for what we mean by exclusive. 
[42] This expression is for a full azimuthal coverage for the lepton detection. For a finite azimuthal 

coverage A</>, it should be multiplied by A4>/2tt. 
[43] The normalization h.P = 1 is in fact irrelevant. The choice h.P=l is a matter of convenience. 
[44] We remind that U u d simply projects on the u, d part of the flavour multiplet 
[45] In the case of pion twist-2 production T also depends on x and on the Sudakov vector n but 
we omit them as they are not important for what follows. 



[46] By definition the argument of spinor is always on the mass shell so we do not specify it. 

[47] Note that the same argument has already been used to neglect the terms shown on|Hfc in the 

saturation of the commutator. Those terms would effectively contribute to the pion poles of 

Fig. GJi and Fig. Ub 

[48] That would not be true for instance in pion photo-production because at threshold t is of order 

2 

[49] Note that there are four N — > A helicity amplitudes for the vector operator. However, gauge 
invariance for the electromagnetic current operator leads to only three electromagnetic form 
factors. Hence there are strictly speaking four GPDs for the vector N — > A transition, of 
which one has a vanishing first moment. In view of our strategy to keep only the dominant 
transitions, we will neglect this GPD which has a vanishing first moment in the following. 

[50] Note that the factor is conventionally chosen in Eq. l(%2*|l such that the Adler form factors 
Cf-(A 2 ) correspond with the p — > A + transition. 

[51] Note that all other (sub-dominant) GPDs vanish at leading order in the l/N c expansion. 

[52] Note that in the large iV c limit, the isovector combination H u — H d is suppressed, therefore 
one could as well give as estimate 2G* M {Q) ~ /j, v ~ 5.43, which is accurate at the 10 % 
level. 

[53] To simplify the notations, we note ABH + ADVCS for the cross sections of the coherent sum 
of both processes. 



